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We investigate the breakup of the one-neutron halo nuclei 11 Be and 
19 C within a dynamical model of the continuum excitation of the projectile. 
■ The time evolution of the projectile in coordinate space is described by solv- 

ing the three-dimensional time dependent Schrodinger equation, treating the 
projectile-target (both Coulomb and nuclear) interaction as a time dependent 
external perturbation. The pure Coulomb breakup dominates the relative en- 
ergy spectra of the fragments in the peak region, while the nuclear breakup is 
important at higher relative energies. The coherent sum of the two contribu- 
tions provides a good overall description of the experimental spectra. Cross 
sections of the first order perturbation theory are derived as a limit of our dy- 
(Nj | namical model. The dynamical effects are found to be of the order of 10-15% 

for the beam energies in the range of 60 - 80 MeV/nucleon. A comparison 
of our results with those of a post form distorted wave Born approximation 
shows that the magnitudes of the higher order effects are dependent on the 
! theoretical model. 
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I. INTRODUCTION 



11 Be and 19 C are examples of one- neutron halo nuclei, where the loosely bound valence 
neutron has a large spatial extension with respect to the respective cores (see, e.g., 
for a recent review). Breakup reactions, in which the valence neutron is removed from the 
projectile in its interaction with a target nucleus, have played a very useful role in probing 
the structure of such nuclei @-10|. Strongly forward peaked angular distributions for the 



neutron p| , [12|j and the narrow widths of the parallel momentum distributions [|l3|-|T9f of the 
core fragments are some of the characteristic features of the breakup reactions induced by 
these nuclei, which provide a clean confirmation of their halo structure. In the Serber types 
of models |i~8l , |20| , the breakup cross section is directly related to the momentum space wave 
function of the projectile ground state. 

The Coulomb breakup is a significant reaction channel in the scattering of halo nuclei 
from a heavy target nucleus (see, e.g., pifl). It provides a convenient way to put constraints 
on the electric dipole response of these nuclei f22| , |2"3|| . The Coulomb breakup of weakly bound 
nuclei can also be used in determining the cross sections of the astrophysically interesting 
radiative capture reactions |[24|| . 

The breakup of the halo nuclei have been investigated theoretically by several authors 
using a number of different approaches (see, e.g., [ID] for an extensive list of references). 



Some of these models [|T^.(||.[26|] use semiclassical geometrical concepts (mostly Serber type) 
to calculate the breakup cross sections. A direct breakup model (DBM) (which reduces to 
the Serber model in a particular limit p7|) has been formulated within the framework of 
the post form distorted wave Born approximation (DWBA) [^8[|l0j. However so far only the 
Coulomb breakup of the halo nuclei has been investigated within this theory. The nuclear 
breakup, for these cases, has been studied mostly within the semiclassical J29| and eikonal 
models 0,0. 

Breakup reactions of halo nuclei can also be described as the inelastic excitation of the 
projectile from its ground state to the continuum [[32|,[3t|. In an extensively used theoretical 
approach, the corresponding T- matrix is written in terms of the prior form DWBA [28]. For 



the pure Coulomb breakup, the semiclassical approximation of this theory is the first order 
perturbative Alder- Winther theory of Coulomb excitation |34j] . It has recently been used to 
analyze the data on the breakup reactions induced by n Be M and 19 C ||. However, higher 
order excitation effects |35|-p9|] may be substantial in the breakup of these nuclei. It has been 
shown [10,40| that for the Coulomb breakup of 19 C, the results of the first order semiclassical 
Coulomb excitation theory differ strongly from those of the DBM which includes the higher 
order effects. 

The general methods in the semiclassical description of the excitation process which 

41~1 , explicit inclusion of 



include higher order effects, are the coupled channel approaches 
successive higher order terms [I2~|, and the direct numerical integration of the time dependent 
Schrodinger equation p6|- p8| , ^3| . The last method, in which all the higher order effects in the 
relative motion of the breakup fragments are included, provides a fully dynamical calculation 
of the projectile excitation caused by both the Coulomb and the nuclear interactions between 
the projectile and the target. 

The time dependent Schrodinger equation method has been used earlier p6| , |37| to inves- 
tigate the Coulomb breakup of n Li and 11 Be. However, these calculations have employed 
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a straight line trajectory for the projectile motion (an approximation valid at higher beam 
energies) and have ignored the spins of the particles in the nuclear models. Furthermore, 
they are partly perturbative in the sense that the first order perturbation theory has been 
used therein to calculate the energy distribution of the breakup cross section at larger impact 
parameters. This procedure may lead |38| to a significantly larger cross section as compared 
to that obtained in models where the partly perturbative approximation is not used. 

In this paper, we present calculations for both the Coulomb and the nuclear breakup of 
the one-neutron halo nuclei n Be and 19 C within the time dependent Schrodinger equation 
method (this will be referred as the dynamical calculation in the rest of this paper). The spins 
of the particles are included explicitly in the nuclear model. The Hamiltonian describing 
the internal motion of the projectile contains a spin dependent interaction. No perturbative 
approximation has been made in the calculation of the breakup cross sections. We also 
calculate first order perturbative results as a limit of our dynamical model so that a consistent 
investigation of the role of the dynamical effect in the breakup cross sections is possible. We 
compare the results of our calculations with those of the post form finite range DWBA theory 
to study the magnitudes of the higher order effects calculated in different approaches. 

Our manuscript is organized in the following way. In section II, we give the details of 
our formalism. The comparison of our calculations with the experimental data and the 
discussion of the results are presented in section III. The summary and conclusions of our 
paper are given in section IV. 



II. FORMALISM 

In this section, we describe a general method for solving the time dependent Schrodinger 
equation, which can be used in problems depending upon a single three-dimensional dynam- 
ical variable. The potentials involved are assumed to be local. The time development of the 
wave function of the relative motion between the neutron and the core nucleus, [^(r, t)], is 
given by the time dependent Schrodinger equation 

ih^(r, t) = [F (r) + V(r, t)} *(r, t), (1) 

where H (r) is the internal Hamiltonian describing the relative motion between the valence 
neutron (n) and the core nucleus (6). It is defined as 

H (r,t)=-^-A + V (r), (2) 

where /if> n is the reduced mass of the n + b system and r = r n — r&. The potential Vb(r) could 
contain a central and a spin-orbit term. In Eq. (p]), V(r, t) is the time dependent external 
field exerted by the target on the projectile. In the present application, it is the Coulomb 
and nuclear interaction between the target (A) and the projectile (a). We write 

*M) = E^^M c (f), 0) 

c 

with 
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yf c M c (*) = E (tcm e sm a \J e M c )Ye cmt (t)x. 



(4) 



mim s 



where c stands for the individual channels characterized by the quantum numbers J c , M c , 
and £ c . The orbital angular momentum (£ c ) is coupled to the neutron spin (s = |) to give 
the total angular momentum (J c ) under the assumption of zero spin for the core nucleus. 
Substituting Eqs. (Q) and @ into Eq. (||), we can write a set of coupled equations for the 
radial wave function 



th—ip c/ {r,t) 



(5) 



where 



h cc '(r,t) 



h 2 



dr 2 



+ V n c (r) 



5 CC , + V cc \r,t). 



(6) 



The solution of the radial equation (|5|) can be written as 

4> c {r,t) = Mr,U) + 4 [ dt ' Y, h cc{r,t')t(j c ,{r,t'), 



(7) 



with the starting time U. Alternatively, the wave function at time t can be calculated from 
the initial wave function by applying the unitary time evolution operator (in the obvious 
matrix notation) 



il>(t)=U{t i ,t)i/>(t i ). 
For a small time step At, the time evolution operator can be approximated as 



XJ(t,t + At) 



1 + ¥nHt) 



(8) 



(9) 



where the elements of the matrix h are given by h cc >. Eqs. (^|) and (||) are the starting points 
for the numerical solution of the problem, which has been done by following the method 
described in Ref. [43]. This approach has the virtue that a first order calculation can also 
be performed within the same program by keeping the coupling only between the initial 
channel (c' = Cj n ) and all the possible final channels (c) in the perturbation potential, i.e., 
by replacing V cc by V cc S c r Cm in Eq. ([]). Then the matrix h is no longer hermitian and the 
time evolution operator U becomes non-unitary. 

Unlike the previous calculations P^55[, we do not assume the projectile to move on a 
straight line trajectory. In our case, the cm. of the projectile is supposed to follow a hyper- 
bolic trajectory with respect to the target during the scattering process. For application to 
the breakup of the one-neutron halo nuclei, the Coulomb part of the external perturbation 
[Vc(r, t)} is given by 



V c (r,t) 



ZAZbe 2 ZAZbe 2 
|r 6 -R A (t)| ~ [RaWI 



(10) 
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where Z A and Z b are the charge numbers of the target and the core nuclei, respectively. 
Rvi(t) is the coordinate of the target in the projectile centre of mass (cm.) frame, and 

r, are the position vectors of the core nucleus and the 



— m^_ r anc i r 



m n +m b 



neutron, respectively. The nuclear part of the external perturbation [Vjv(r, £)] is the sum of 
the neutron-target and core-target optical model potentials, respectively. It is given by 



V N (r,t) = V n f(\r n - R A (t)\,R Vn ,a Vn ) +iW n f(\r n - R A (t)\,R Wn ,a Wn ) 
+ V b f(\r b - R A (t)\,R Vb ,a Vb ) + iW b f(\r b - R A (t)\, R Wb ,a Wb ), 



where 



f(x,R,a) 



1 + exp 



x — R 



(11) 



(12) 



with the diffuseness parameter a and radius R. The depths of the real and the imaginary 
parts of the neutron-target and the core-target optical potentials are denoted by V n , W n , 
and V b , W b , respectively. The corresponding radii and the diffuseness parameters are given 
by Ry n , Rwm and ay n , o-Wn, and Ryb, Rwb and a Vb , awb, respectively. 

Using a multipole expansion of the perturbation potentials, we can write 



V-\r 1 t)=Y.C%[v^{r,t)+v x N ^t) 



with the coefficients 



dQ 



y 



JcMc 



The time dependent radial potentials are given by 



— 77V. 



^KZ A Z b & 

2A + 1 \m n + m b 



K A 



(*) 



(13) 



(14) 



(15) 



and 



W> t) = V n f x ^(r n , R Vn , avn, t) + iW n f x ^(r n , R Wn , a Wn , t) 
+ V b f Xf *(r b , R vb , a Vb , t) + iW b f x »(r b , R Wb ,awb,t), 



where 



f x »(r, R, a,t) = J dQ /(|r - R A (t)\,R, a) Y*^ 



r . 



(16) 



(17) 



The triple differential cross section for the breakup of the projectile a into b and n is 
given by 



d a 



a 



da 



R 



dE bn dVt bn dVt bn - A dQ 



bn-A 



-Pif(kbn) Pf(hn), 



(18) 



where Q bn -A are the angles associated with the relative motion of the cm. of the b + n 
system with respect to the target nucleus A. The first term on the right hand side of Eq. 
( f[8|) is the Rutherford cross section for the projectile-target scattering. The last term, 
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Pbnhn (lg) 

Pf (2tt)^ 2 ' 1 j 

is the density of final states. The excitation probability, (which is a function the relative 
momentum Kk bn between the two fragments), is given by 

P tf (M= lim — i- £ kCHk^r)!^^^/)}! 2 , (20) 

where $^(k bn , r) is a complete scattering solution for the relative motion of the two frag- 
ments with the in-going wave boundary condition. It satisfies the Schrodinger equation (|]) 
for a vanishing perturbation potential V(r,t). The wave function ^^(r, tf) is the solution 
of the time-dependent Schrodinger equation with the boundary condition that as U — > — oo, 
it goes to the unperturbed wave function $^ of the ground state of the projectile. Since we 
are interested only in the breakup contributions, it is advantageous for numerical reasons to 
replace the full wave function ^^(r, £), in Eq. (|2~0"|), by the corresponding continuum wave 
function 

*S?(r,t)=* JMi (r l t)-X;<^(r,t)|*«i(r l t)>$y(r,t), (21) 

v 

where the sum runs over all bound states of the system. 

In the actual numerical realization, the radial wave functions ipc{ r ,t) are discretized on 
a mesh with points x n = nAx where Ax = 0.0025 and n — 0, . . . , 400. They are related to 
the points in the radial coordinate r by the mapping, r n = /2 mcta; [exp(aa; n ) — l]/[exp(a) — 1], 
with R m ax — 900 fm. The parameter a is chosen so that r\ = 0.3 fm. The second order 
derivative, J^, in Eq. (^j) is represented by a finite difference approximation. In the first 
step, the ground state radial wave function is calculated assuming a Woods-Saxon potential, 
Vq{t) = Vof(r,R ,a )5 cc >, for the n — b interaction in the initial channel with a given set 
of quantum numbers, J c , M c , and l c . The time evolution of this initial wave function is 
calculated by the repeated application of the time evolution operator @, using time steps 
At = 1 fm/c, and taking into account partial waves with I = 0, 1, 2, 3 for the n — b relative 
motion. This leads to a set of coupled linear equations for the radial wave functions which 
are solved with the technique as described in Ref. []43[1 . We use a time interval for the 
evolution which is symmetric to the time of closest approach between the target and the 
projectile. Its limits are determined by the condition that the perturbation potential is at 
least 200 times smaller than its maximum value. In order to avoid spurious excitations, the 
time dependent potential was switched on adiabatically. Furthermore, care has been taken 
to ensure that unphysical bound states are not populated during the time evolution. 

TABLE I. Optical potential parameters for the neutron-target and core-target interaction. 



projectile 


V n [fm] 


R Vn [fm] 


avn [fm] 


W n [fm] 


R Wn [fm] 


a Wn [fm] 


Ref. 


n B 


-27.88 


6.93 


0.75 


-14.28 


7.47 


0.58 


a 


19 C 


-29.48 


6.93 


0.75 


-13.18 


7.47 


0.58 


a 


projectile 


V b [fm] 


Rvb [fm] 


a Vb [fm] 


W b [fm] 


Rwb [fm] 


a Wb [fm] 


Ref. 


n Be 


-70.0 


5.45 


1.04 


-58.9 


5.27 


0.887 


m 


19 C 


-200.0 


5.39 


0.90 


-76.2 


6.58 


0.38 


n 
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The final wave function is projected onto scattering states and the excitation probability is 
calculated according to Eq. (^Of), after removing the bound state contributions. For each 
M c substate of the initial state, an independent calculation has been performed. 

III. RESULTS AND DISCUSSIONS 

The parameters of the nuclear optical potentials for the n-target and core-target interac- 
tion [see, Eq. (16)] used in our calculations, are given in Table I. We have adopted a single 
particle potential model to calculate the ground state wave function of the projectile. The 
ground state of n Be was assumed to have a 2s\/2 valence neutron coupled to the + 10 Be 
core with a binding energy of 504 keV. The corresponding single particle wave function was 
constructed by assuming the neutron - 10 Be interaction of the Woods - Saxon type having a 
central and a spin-orbit part. The radius and the diffuseness parameters in both the terms 
were taken to be 2.478 fm and 0.5 fm, respectively. The depth of the central term was 
searched so as to reproduce the ground state binding energy. This 2si/2 wave function has 
an additional node as compared to a simpler zero-range wave function. The strength of 
the spin-orbit term was adjusted by requiring that the first excited (~ ) state in n Be is at 
the experimental excitation energy of 320 keV. In the dynamical calculation, we assume a 
vanishing n-core potential for the higher partial waves so that unphysical resonances can be 
avoided. For 19 C, the ground state wave function was obtained with the similar procedure 
by assuming a configuration in which a 2sx/2 neutron is coupled to the + 18 C core. The 
radius and the diffuseness parameters associated with the Woods-Saxon interaction were 
taken to be 3.30 fm and 0.65 fm, respectively. In the dynamical calculations, we took into 
account only the n-core interaction in the s-wave. The spectroscopic factors for the ground 

2.5 I I 
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FIG. 1. The differential cross section as a function of the relative energy of the fragments 
(neutron and 10 Be) emitted in the 11 Be induced breakup reaction on a 208 Pb target at the beam 
energy of 72 MeV/nucleon. The long-dashed and short-dashed lines represent the pure Coulomb 
and pure nuclear breakup contributions, respectively. Their coherent sum is depicted by the solid 
line. The experimental data are taken from ||. 
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state were taken to be 1 in all the cases. In the multipole expansion of the perturbation 
potential [Eq. (|13|)1, we included contributions up to A = 2. 

In Fig. 1, we compare the results of our calculations with the experimental data (taken 
from H) for the relative energy spectrum of the fragments emitted in the breakup reaction 
of n Be on a 208 Pb target at the beam energy of 72 MeV/nucleon. In these calculations, the 
integration over the theta {6bn-A) angles of the projectile cm. was done in the range of 0°— 3°, 
which corresponds to a minimum impact parameter of about 12 fm. The long-dashed and 
short-dashed lines represent the results of the dynamical calculations for the pure Coulomb 
and the pure nuclear breakup, respectively, while their coherent sum is represented by the 
solid line. The Coulomb cross sections consist mostly of the dipole term as the contributions 
of the quadrupole mode are negligible. 

We note that the pure Coulomb contributions dominate the cross sections around the 
peak value while the nuclear breakup is important at the larger relative energies. This can 
be understood from the energy dependence of the two contributions. The nuclear breakup 
occurs when the projectile and the target nuclei are close to each other. Its magnitude which 
is determined mostly by the geometrical conditions, has a weak dependence on the relative 
energy of the outgoing fragments beyond a certain minimum value. Contrary to this, the 
Coulomb breakup contribution has a long range and it shows a strong energy dependence. 
The number of virtual photons increases for small excitation energies. At the same time, 
a much larger range in the impact parameter implies a substantial breakup probability. 
The results shown in this figure are in agreement with those of Ref . ||4| . The domination 



of the nuclear breakup may explain the failure of the pure Coulomb finite range DWBA 
calculations |10[ in explaining the data at larger relative energies. 

The coherent sum of the Coulomb and the nuclear contributions provides a good overall 
description of the experimental data. Although, the pure nuclear breakup contributions 
below 0.5 MeV are substantial, yet their interference with the Coulomb part does not change 
appreciably the shape of the total cross section in the peak region. It may be noted that 
the absolute magnitude of the cross section near the maximum is somewhat underestimated 
by our calculations, but the position of the peak is well reproduced. This is in agreement 
with the results of the dynamical calculations reported in . However, in the semiclassical 
coupled-channel calculations |Q for this reaction, the peak position of the calculated cross 
sections are shifted towards larger energies as compared to that of the data. 

It may be remarked that use of a some what smaller minimum impact parameter in the 
angular integration over the n Be cm. would increase the cross sections which may reduce 
the difference between the experiment and the theory. However, this would mainly affect 
the nuclear contribution since it is more sensitive to the smaller projectile-target distances. 
The Coulomb contributions arise mostly from the much larger impact parameters as long 
as the excitation energy is not too large. At very small impact parameters the absorption 
due to the imaginary parts of the optical potentials will result in an effective cutoff. Since 
the optical potentials are not sufficiently well known, the choice of the maximum scattering 
angle introduces an uncertainty but our choice for this is a reasonable one. 

In Fig. 2, we show the results of our dynamical calculations for the relative energy 
spectrum of the fragments (neutron and 18 C) emitted in the breakup of 19 C on a 208 Pb 
target at the beam energy of 67 MeV/nucleon. The integrations over Qbn-A were done in 
the range of 0° — 3°. Since the binding energy (B n _ core ) of the valence neutron - core system 



S 



in the ground state of 19 C is still an unsettled issue [j^f46[ -|49|j , we present in this figure the 
results of calculations performed with two values [530 keV (part a) and 650 keV (part b)] 
of B n _ core . This is for the first time that the dynamical calculations (including the nuclear 
breakup) have been performed for this case. The experimental data are taken from ||. We 
see that the Coulomb breakup dominates the cross sections in the peak region while the 
nuclear breakup is important at the larger relative energies in this case too. 
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E rel [MeV] 

19 C + 208 Pb-> 18 C + n + 208 Pb (b) : 

E beam = 67 MeV/nucleon 

Coul. dynamical j 

Coul. + Nucl. dynamical 

Nucl. dynamical 

"'"O.O 0.5 1.0 1.5 2.0 2.5 3.0 

E rel [MeV] 

FIG. 2. The relative energy spectrum of the fragments (neutron and 18 C) in the breakup of 19 C 
on a 208 Pb target at the beam energy of 67 MeV/nucleon. Part (a) shows the results obtained with 
a binding energy of 530 keV for the neutron-core configuration in the ground state of 19 C while 
part (b) is the result with a value of 650 keV for the same. The dashed and dotted lines represent 
the pure Coulomb and pure nuclear breakup contributions, respectively, while their coherent sum 
is shown by the solid line. The experimental data are taken from p| . 

The dynamical calculations carried out with the binding energy of 530 keV [part (a)] 
overestimate the cross sections in the peak region by about 35 - 40%, while those done 
with 650 keV [part (b)] are in good agreement with the data. Our calculations, therefore, 
seem to support the latter value for B n _ core which is within the error range of the value 
[(530 ± 130) keV] reported in |§. The expected shift in the peak position in part (b) is 
within statistical error of the data. We would like to recall, however, that we have used a 
spectroscopic factor of 1 for the neutron-core configuration for the ground state of 19 C. 

An alternative scenario has been presented in |J , where the same data have been analyzed 
within the first order semiclassical perturbation theory of the Coulomb excitation. With a 
binding energy of 530 keV and the same neutron-core configuration for the ground state of 
19 C, these calculations overestimate the data in the peak region by about 40-50% (see also 
p0|). Instead of using a larger value of B n __ core in order to fit the data, these authors try to 
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extract a spectroscopic factor (SF) for this configuration by comparing their pure Coulomb 
dissociation calculations with the data (which are corrected for the nuclear breakup effects 
in an approximate way). The value of SF (0.67) obtained in this way is close to the shell 
model results of Ref . [|5| . 

However, a note of caution must be added about the significance to this result. The 
spectroscopic factors reported in ||5|] correspond to various specific states of the 18 C core. 
The data of Ref. || are, however, inclusive in the sense that the specific core states of 18 C are 
not observed there. Therefore, measurements of the relative energy spectra in experiments 
of the type reported in pEfl , where specific core states of 18 C are identified by tagging to 
the decay photons, are required for a meaningful comparison of the SF extracted from such 
studies with those of the shell model. Moreover, the procedure adopted in || to correct the 
data for the nuclear breakup effects may not be valid due to the long range of the nuclear 
interaction in case of the halo nuclei [[4ij| . Obviously, the value of the SF extracted from the 
breakup studies depends critically on the theory used to calculate the corresponding cross 
sections. Therefore, more experimental data and their theoretical analysis within different 
models of the breakup reactions are required for arriving at a more definite conclusion in 
this regard. 
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FIG. 3. Dynamical model results for the angular distribution of the center of mass of n + 18 C 
system in the breakup of 19 C on a Pb target at the beam energy of 67 MeV/nucleon for two values 
[530 keV (part a) and 650 keV (part b)] of the 19 C ground state. The dashed and dotted lines 
represent the pure Coulomb and pure nuclear contributions. Their coherent sum is represented by 
the solid line. 



Another argument put forward in || in favor of their procedure is that with the values 
0.530 MeV and 0.67 for B n _ core and SF, respectively, the angular distribution of the n + 18 C 
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cm. measured in the same reaction is well reproduced. However, this conclusion assumes 
that the shape of this angular distribution is not affected by the nuclear breakup effects 
below the grazing angle (~ 2.7°). To study the role of the nuclear breakup for this data, we 
show in Fig. 3 the angular distribution of the cm. of the n + 18 C system for the values of 
B n - core of 530 keV (upper part) and 650 keV (lower part). The integrations over the relative 
energy is performed in the range of 0.0-0.5 MeV (the same as done in Ref. 0). It is clear 
from this figure that the nuclear breakup effects start becoming important already from 
cm. angles of 1.5°, and the total cross section (coherent sum of the Coulomb and nuclear 
contributions) differs from the pure Coulomb one even below the grazing angle. Beyond this 
angle, the absorptive part of the optical potentials reduces the cross section. In this figure 
we have not shown the comparison of our calculations with the experimental data as it would 
require folding the calculations with the experimental angular resolution II. However, the 



comparison of the unfolded and folded results as shown in |4(| suggests that the quality of 
agreement between calculations done with both the options and the data would be similar. 

The magnitude of higher order effects in the Coulomb breakup reactions of n Be and 
19 C is a subject of current interest. In an earlier calculation |[36| , higher order effects were 
found to be rather small for the reaction studied in Fig. 1. However, comparing the result of 
the adiabatic model of Coulomb breakup reactions with that of the first order semiclassical 



perturbation theory of the Coulomb excitation, it has been concluded in [|39] that the higher 
order effects are substantial for the reaction investigated in Fig. 2, which 



9_ 208-. 18- 208i-ii_ 

C + Pb -> C + n + Pb 



E beam = 67 MeV/nucleon 



Coul. first order, 530 keV 
Coul. dynamical, 530 keV 

Coul. first order, 650 keV 

Coul. dynamical, 650 keV 




[MeV] 



FIG. 4. Comparison of the dynamical model and the first order perturbation theory results for 
the pure Coulomb breakup contribution to the relative energy spectrum of the fragments (neutron 
and 18 C) emitted in the breakup of 19 C on a 208 Pb target at the beam energy of 67 MeV/nucleon. 
The solid and dotted lines represent the results of the first order and the dynamical model calcu- 
lations, respectively, obtained with the value of B n _ core equal to 530 keV while the dashed and 
dashed-dotted lines represent the same for the B„,_ core value of 650 keV. The first order perturba- 
tion theory results have been obtained from the dynamical model in a particular limit as discussed 
in section II. 



11 



would have a considerable influence on the extracted spectroscopic factor. In contrast, these 
effects were found to be rather small 



in 



50 



where the first order and the higher order 
terms were calculated within the same model. As discussed in section II, the first order 
perturbation theory results can be obtained from our dynamical model in a particular limit. 
In Fig. 4, the first order results obtained in this way are compared with those of the full 
dynamical model for the pure Coulomb breakup contribution to the same reaction as in Fig. 
2. It can be noted that the higher order effects are rather small. They reduce the peak 
cross sections of the first order theory by about 10% for both values of B n _ core but leave the 
shapes of the spectra largely unaffected. 

In Fig. 5, we compare the pure Coulomb breakup contributions (calculated within the 
dynamical model and the finite range DWBA [|Kj) to the reactions studied in Figs. 1 and 2. 
In the finite range DWBA theory the higher order effects in the target-fragment interaction 
are automatically incorporated. The interesting aspect of this comparison is that while for 
the 11 Be case the finite range DWBA results are slightly larger than those of the dynamical 
semiclassical calculation, the former is smaller than the latter for 19 C induced reaction for 
both the values of the binding energy. This clearly shows that higher order calculations for 
the Coulomb breakup performed within different theories could be different from each other. 
The reason for this difference is the fact that fully quantal higher order approaches 




ir, 208,-.. 10,-. 208 — . 

Be + Pb -> Be + n + Pb 
E beam = 72 MeV/nucleon 



Coul. dynamical 
FRDWBA 



[MeV] 



19„ 208,,. 18„ 208-. 

C + Pb -> C + n + Pb 
E beam = 67 MeV/nucleon 




(b) 



Coul. dynamical, 530 keV 

Coul. dynamical, 650 keV 
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FRDWBA, 650 keV 
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FIG. 5. (a) Comparison of the pure Coulomb dynamical model (solid line) and the finite range 
DWBA (dotted line) contributions to the relative energy spectrum of the fragments (neutron and 
the core) emitted in the pure Coulomb breakup of 11 Be on a 208 Pb target at the beam energy of 
72 MeV/nucleon. (b) The same as in part (a) for the breakup of 19 C on the same target at the 
beam energy of 67 MeV/nucleon. In this case the full and dashed lines represent the results of the 
the dynamical model and finite range DWBA, respectively, corresponding the 19 C binding energy 
of 530 keV while dashed and dashed-dotted lines show the same for the binding energy of 650 keV. 
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take into account effects which are beyond the semiclassical dynamical model calculations. 
Therefore, the results of the semiclassical first order perturbation theory will differ from those 
of different higher order models in different ways. Thus, the conclusion about the role of the 
higher order effects will differ if the first order and higher order calculations performed within 
two different theories are compared with each other. It may also be remarked here that the 
spectroscopic factors extracted from the comparison of the Coulomb breakup calculations 
of different models with the corresponding data could be quite different from each other, as 
has already been pointed out in j4(J. 

IV. SUMMARY AND CONCLUSIONS 

In this paper, we have studied the Coulomb and the nuclear breakup of the one-neutron 
halo nuclei, n Be and 19 C, on a 208 Pb target within a semiclassical fully dynamical model of 
the projectile excitation process. The time evolution of the projectile system is described by 
numerically solving the time dependent Schrodinger equation in three dimensions. Through 
this nonperturbative method higher order effects are fully taken into account in the breakup 
process. Unlike previous such calculations, we assume the projectile to move along a hyper- 
bolic trajectory instead of a straight-line one. A simple single particle potential model was 
used to calculate the ground state wave function of the projectile nuclei. 

With a configuration for the n Be ground state in which a 2s\/2 neutron is coupled to 
the ground state of the 10 Be core with a binding energy of 0.504 MeV and a spectroscopic 
factor of unity, the experimental relative energy distributions of the fragments (neutron and 
10 Be) emitted in the breakup reaction of n Be on a lead target at the beam energy of 72 
MeV/nucleon are described rather well by our model. The Coulomb breakup dominates 
these cross sections around the peak region while the nuclear breakup is important at higher 
relative energies. This provides a natural explanation for the failure of the pure Coulomb 
breakup calculations in describing these data at the larger relative energies. 

For the 19 C case, the comparison of the dynamical calculations with the data for the 
relative energy spectrum of the fragments suggests that the ground state configuration of 
this nucleus is consistent with a 2si/ 2 neutron coupled to the ground state of the 18 C core 
with a binding energy of 0.650 MeV and a spectroscopic factor of 1. Alternatively, it could 
also have a binding energy of 0.530 MeV but a spectroscopic factor of about 0.6 — 0.7. 
The current data || on the relative energy spectrum are not sufficient to rule out either of 
these possibilities. Further experimental studies in which the relative energy spectrum is 
measured by tagging the specific core states and calculations of the Coulomb and nuclear 
breakup effects within different models (e.g., finite range DWBA) are clearly needed to settle 
this issue. 

By comparing the dynamical model results for the pure Coulomb breakup with the first 
order calculations (which are deduced from the same model in a particular limit), we find 
that higher order effects in these reactions are generally small. They reduce the magnitudes 
of the first order relative energy distribution of the fragments by about 10-15% in the peak 
region. Their shapes, however, remain almost unaltered by the higher order effects. However, 
a comparison of the higher order and first order calculations of two different models may lead 
to large differences in the two results. This is due to the fact that higher order calculations 
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performed within different theories include these effects in different approximations and 
could differ quite a bit from each other. 
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